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, Abstract. We investigate a semilinear modification for the wave map prob- 

l^^i lem proposed by Adkins and Nappi [I], and prove that in the equivariant case 

the solution remain continuous at the first possible singularity. This is usually 
one of the steps in proving existence of global smooth solutions for certain 
equations. 



1. Introduction 



< 

''pi I Let (j) '■ IR"+^ — > M be a map from the n + I dimensional spacetime, with 

C^ ' Lorentzian metric g of signature {l,n), to a Riemannian manifold {M,h). The 

action of the wave map equation, or the nonlinear a model in physics terminology, 

is 

> . (1) ^ ^ ^ / -^""^ ^^'^^ ^-^' ^^^" ("^^ "^9 

Q\ , The initial value problem for the Euler-Lagrange equations associated with (|T]) 

^D ■ has been intensely studied, especially the issues of global existence and regularity for 

its solutions. We mention here the pioneering works of Christodoulou - Tahvildar- 

lO ■ Zadeh [5 , Grillakis [8 , Shatah - Tahvildar-Zadeh [l^, and Struwe [16]. 

^P , The particular case when M = S'^ and n = 3 is of special interest in high energy 

physics. The nuclei of atoms are held together by forces mediated by the pi mesons. 
These are a set of three particles whose masses are small compared to the nuclei 
themselves, so to a first approximation they can be considered to be massless, 
i.e., travelling at the speed of light. If interactions among them are ignored, the 
pi mesons are described by a field (j> : R^'^^ — >■ R^ satisfying the wave equation. 
5^ I Interactions would add nonlinearities. A remarkable fact of physics is that the 

interactions among the pi mesons are described, to a good approximation (Gell- 
Mann - Levy [7] , Gursey [5] , [TU] , and Lee JT] ) , by considering the target manifold 
to be the sphere S'^ and replacing the wave equation by the corresponding wave 
map equation. Physicists call this the nonlinear a model for historical reasons. 

In order for the energy to be finite, the gradient of the field must vanish at 
infinity: it tends to the same value in all directions of spatial infinity. Thus, each 
instantaneous pi meson configuration corresponds to a map </) : S"^ — > S"^ obtained 
by identifying the points at spatial infinity. A continuous map of the sphere to itself 
has an integer associated with it, the winding number or element of the homotopy 
group 7r3(S'^). In terms of the original field, we can write this as 
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(2) Q = cJ e,,kda^'db^d,(j)^e"-'^ dx 

where e is the Levi-Civita symbol and c is a normahzing constant. Small perturba- 
tions from the constant configurations would have (5 = 0. If time evolutions were 
continuous, Q would be a conserved quantity, a 'topological charge'. 

It was suggested by Skyrme [ini in the 1960s that this topological charge Q is 
just the total number of neutrons and protons in a nucleus: the baryon number. 
Thus, a proton (hydrogen nucleus) would have Q = 1, a helium nucleus would 
have either Q = 3 or 4 (depending on the isotope), and so on. This has been 
confirmed by connecting with Quantum chromodynamics (Balachandran - Nair 
- Rajeev - Stern [2], [3], Witten [18], [H]), the fundamental theory of nuclear 
interactions. Numerical calculations (Battye - Sutcliffe [1]) also support this idea. 
Some essential theoretical puzzles had to be resolved before this rather strange idea 
of Skyrme could be established. These were resolved in the mid-80's; e.g., by using 
topological notions ('anomalies') ([E], 0) [2], [E], [19]). 

One of these puzzles is that actually the wave map does not have continuous time 
evolution: it is supercritical. Shatah |13| has exhibited an example of a solution 
with smooth initial conditions that breaks down in finite time. Physically, this is 
because the forces among the pi mesons are mainly attractive; so a configuration 
of winding number one would shrink to a point, emitting its energy as radiation 
carried away to infinity. At the singularity, Q would jump to zero. Since baryon 
number is strictly conserved by nuclear forces, this cannot be right. 

Skyrme suggested a modification of the action of the wave map that could sta- 
bilize configurations with Q ^ 0. The corresponding equations are quite intricate, 
being quasilinear. What we study here is a semilinear modification of the above 
action S", previously introduced by Adkins - Nappi jl], which is also expected, 
for physical reasons, to have regular solution^. The idea is to add a short range 
repulsion among the pi mesons, created by their interactions with an omega me- 
son. Geometrically, it describes the nonlinear coupling of the wave map with a 
gauge field, the source (charge density) of the gauge field being the density of the 
topological charge Q. 

In detail, the action of the theory we investigate is 

(3) S ^ S + ^J F^'-'F^, dg - J A^r dx 

where the omega meson is represented by the gauge potential A = Af^dx^^, the 
2-form Ffii, = d^Ai, — duA^ is its associated electromagnetic field, and 

is the fiux or the baryonic currenlo. 

Here we study time dependent equivariant maps associated to S with winding 
number 1, i.e., 

g = diag(-l, 1,1,1) A/ = §3 
0(i, r, ^, 9) = {u{t, r), ^, 0) A{t, r, V^ 9) = {V{t, r), 0, 0, 0) 



We consider only the case when the masses of the mesons are set to zero, which is sufficient 
to understand short distance singularities. 
^Note that Q = J f dx. 
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with boundary conditions u{t,0) = 0, u{t, oo) — tt, and V{t,oo) = 0. 
Routine computations lead to 



2n 



9, , ,, „ . 9 a^(w — sinucosw)^"' 

H i-uf + m;) + 2 sm^ u H ^ ^ 



drdt- 



rVr 



a{u — sin u cos u)"' 



dr dt 



where a is an appropriate constant. Ehminating V using its variational equation 
and scaling the coordinates by a factor of a, we obtain the main equation satisfied 
by u: 

2 1 1 

(4) utt — Urr Ur H — 9 sin 2u H — jiu — sinucosu)(l — cos2u) = 

Based on the physical intuition detailed above, we conjecture that smooth 
finite energy initial data for this equation will evolve into globally regular 
solutions. This is supported also by the existence of static numerical solutions of 
winding number 1 for ([3]), which were found in jT]. Regularity at the origin forces 
the initial data to vanish at r = 0, which is consistent with the previously imposed 
boundary condition, u{t, 0) = 0. 

One needs to compare ([4]) with the corresponding wave map equation 

2 1 

(5) Utt — Urr Ur H — j sin 2u = 

for which Shatah's solution (later found by Turok and Spergel [JT in close form) 

r 



u{t, r) — 2 arctan ■ 

To — t 

provides an example of smooth initial data that develop singularities in finite time. 
The goal of this article is to take the first step in proving our conjecture, which 
is to show that a smooth solution for ^ remains continuous at the first possible 
singularity. Using the fact that Q is semilinear and it is invariant under trans- 
lations, we can assume, without any loss of generality, that our solution starts at 
time t = —I and breaks down at the origin (0, 0). Our main result is 

Theorem 1.1. The solution u for Q is continuous at the origin and 

u{t,r)->0 as (i,r)^(0,0) 



The proof follows the lines of the one for 2 + 1 dimensional equivariant wave 
maps (e.g., [H]); however, because we are in 3 + 1 dimensions, we lose a critical 
estimate which is in turn bypassed by a new argument, using the sign of the last 
term in ^■. 

(6) u ■ —r(u — sin u cos u) (1 — cos 2u) > 

In all of these problems (i.e., wave maps and Skyrme model) the initial configu- 
ration shrinks as time evolves, causing energy and winding number to accumulate 
at the origin. For the 2 + 1 dimensional wave map equation, enough energy is ra- 
diated away so that the energy density at the origin is finite. In 3 + 1 dimensions, 
the energy density diverges at the origin although the total energy tends to zero. 
For the Adkins - Nappi version of the Skyrme model, there is a repulsive part in 
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the energy density. Our new argument shows that the energy corresponding to this 
part does not concentrate, ahowing us to prove the continuity of the field. 

2. Main argument 

We will be working with backward truncated cones, their mantels and bases, 
denoted as 

i^l = {{t, r)\T<t<S,0<r< \t\} 

C| = {{t,r)\T<t<S,r=\t\} 

Bt = {(i,r)|t = r, 0<r<|i|} 

where —1<T<S<Q. For narrower bases, we use the notation 

Bt{\) = {(t,r)| t^T,-Xt<r< -t} 

where < A < 1. 

Multiplying (U) by r^ut, r^Ur, r^Ur, r^u, respectively r'^tut, we obtain the fol- 
lowing differential identities: 

Lemma 2.1. A classical solution for ^ satisfies: 



(7) dt ( -;r{ut + Uj.) -f sin u + -^ri^{u — sinucosu) ) — dr{r UtUr) — 



(8) 



^2 2^ 

dt{r^UtUr)—dr ( — (wj + Ur) ^ sin u — 77-2(^ ^ sinucosu)^ 



2 ^ ' " 2r2 

1 



r{Uj, — Mt ) T ("" ~ sin u cos u) 



(9) 



dt{r uiUr)~dr — (wj + u^) — r sin u ("" ^ sin u cos u) 

„2 1 

',2 o„,2\ I „;„2„, -'^ /„, „■ „,„ „,^2 



— (u^ — 3uj ) -I- sin u — ;p2 (^ ~ sin u cos u) 



2 ^ '■ ^' 2r2 

„2/„,2 



9t(r uut) — dr{r uur) = r {u^ — u,.)— 

(10) f . 1 

— u sin2u H — k{u — sinucosu)(l — cos2u) 



(11) 



dt i -T- (ui + u^) + t sin u + —^ {u ~ sin u cos u) \ — dr{tr utUr) ~ 

„2 -1 

— —-(Uf + u„] + sni u + —^[u — smucosu) 
2 2r^ 

Remark 2.2. T/ie argument for 2-1-1 dimensional equivariant wave maps relies 
on the counterparts of all five identities. For our analysis we will not use ([8]). 

Next, we define the local energy and the fiux associated to our problem: 

Definition 2.3. The energy of the time slice t — T is defined as 



t / o n-, t . 9 1 

Bt 



E{T) = I ^{u^ +ul) + —sin^u+ —^{u-sinucosu)'^ 
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while the flux between the time slices t = T and t = S is given by 

1 

2^ 



F(T,S) = —= / -{ut - UrY + ^T SVC? u ^ -(w- sinucosu)^ 

V2 7c| 2 r^ 



The classical energy estimate, obtained by integrating ([7]) over K^^ 

E{T)~E{S) = F(T,5) 
implies that the local energy is decreasing and the flux decays to 0: 

(12) E{S) < E{T) for -1<T <S 

(13) lim F(T,0) = 
These two facts allows us to go further and show that: 

Proposition 2.4. The solution u is bounded, with 

(14) ||u|U» < C{E{-1)) 



an- 



d 



(15) lim u{t,-t) =0 

t->o- 



Proof. For 



we have 



I{z) — / (w — sin w cos w) (iti; 
Jo 



2 • 2 

z — sm z 



/(O) = /(z) > 0(z 7^ 0) lim I{z) = oo 

|2:|— ^CXD 



Based on u{t, 0) = we can write 



I(u(t,r)) — / {u — sin u cos u)ur dr 
Jo 



which implies 



T// ^^ ^ I I 7 \' II (u — sinucosu)^\ ^ , ^, , 



iBt J \JBt 

The monotonicity of the energy (fT2|) then settles the first claim. 
Next, denoting v(t) = u{t,—t), we use p^ to infer 

/■° (u — sinwcosw)^ , 
Urn / -^ ^ ^ dt = 

T-^0- Jrp i^ 

which leads to the existence of a sequence t„ — > for which v{tn) -^ 0. For fixed T 
and large n we obtain 

|/(w(T))-/(w(t„))| < / \{v ~ ainvcosv) ■ Vs\ds < F{T,0) 
Jt 

which proves ([T5|). D 

We use these results to reduce the proof of Theorem 11.11 to the one of a local 
energy estimate. 
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Theorem 2.5. // the solution u is smooth on K*_i for all —I < t < and 

16) lim / -^ !- =0 

then u is continuous at the origin and 

u{t,r)-)-0 as (i,r)^(0,0) 



Proof. We argue as in Proposition 12.41 to deduce 

I{u{t,r)) = I{u{t,—t))+ / {u — smucosu){t,r') ■ Ur{t,r') dr' < 

(u — sinucosw)^^ ^ 



<I{u{t,-t))+E{t)2 ■ 

'Bt 

which obviously provides the desired conclusion based on ([T2|) and ([TS]) . D 

Remark 2.6. Theorem \2.5\ is the point where our argument leaves the approach 
for equivariant wave maps. There, one relies on the weaker bound 

I{u{t,r)) < I{u{t,-t)) + E{t) 

and proves nonconcentration of the energy (i.e., E{t) — > as t — > 0/ The crucial 
ingredient for that analysis is that the annular energy doesn't concentrate. 



Exit) = / e^O 

JBtiX) 

where e is the energy density and < A < 1 . This is in turn obtained from 

(17) {{rm)t — (re)r.)^ ^ (e — m){e + m) 

where m is the momentum density. We refer the interested reader to [5] or [12j for 
more details. 

In our case, the corresponding estimate should read 

(18) ((r^m)t - {r'^e)rf < r^(e - m)(e + m) 
where 

e — -iui + ul) H — - aiv? u H -riu — sin u cos u)^ m = UrUt 

2 r^ 2r^ 

The left hand side in (|18|) . which can be obtained through ^, takes the form 



I 1 \ I 1 \ I 1 i\ (w — sinwcosu)^ . 

(r m)t — [r e)r — r[u^. — u^) -\ :: — 2 sm 2u ■ u^ — 

2(u — sinucosw)(l — cos 2m) 
2 "^ 

Precisely the last term above fails to obey the bound in psp. This is the reason why 
we do not use (jH)) in our argument. 
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3. Local energy estimates 

These estimates are deduced by integrating the differential identities (|9|)-(fTT|) on 
the backward cone Kj. and then aUow for T — > 0. First, we notice that for / e L°° 
one obtains 

1 r f ,,„ f f 

r'' 



(19) 



lim 



hm 







which aUows us to ignore certain terms in the analysis 



Lemma 3.1. For u solution of ^, smooth on K*_i for all — 1 < t < 0, the 

following estimates hold: 



(20) lim -i- / 

(21) lim -i- / 

(22) lim -i- / 



KO 



K" 



-(3u^ — u^) + -z-^{u — sin u cos u)^ 



— — ■ / ruj-ut — 

\T\ ' 



Uj, — Uf + u ■ —{u ~ sinMC0sw)(l — cos2u) 



'7T(^f + "r) ^ TT-ri"^ ^ sin M cos w)^ 
2 2r^ 



E{T) = 



Proof. We prove only the first estimate, the other two being treated similarly. As 
mentioned above, we integrate (j9|) on K^ to infer 



K" 



-(3m? — u^) T sin^ u -\ -r(u — sin u cos u)^ 

2 r^ 2r^ 



rUrUt 



+ 






— (uj — Ur)^H — ^ sin u -\ -(u — sinucosu) 



2r4 



Using ([T^ we deduce that 



1 /■ 1 . 2 
hm - — - / —;r sm u = 



(PO)) follows immediately as 



lim - — - 



r ■ I —(ut ~ Ur) H — 7T sm u H j(w — smucosw) 

2 r^ 2r'' 



= 



due to ^ and that r < |T| in C^. 



D 



Finally, combining p0 | -p2 |) and relying on (|6]), we obtain: 

1 



(23) 

Coupling this with 



lim EiT) , , 



rUrUt = 



Bj" 



''^'-wiL 



rUrUt 



B'j' 



-(1 - -^Mut - Urf + ^(1 + ^)(wt + Urf 



1 -.2 1 



Tj sin'^ u H J (w — sin u cos u) 



we deduce the sufficient condition (ITBl) from Theorem 12.51 which finishes the argu- 
ment. 



lim 

T-^0- 



r 1 

\1 \ r^ 



/. 



Bt 1^1 
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It is worth noting that ([23]) yields also 

1 „ 1 

= 

Thus we obtain that the entire energy does not concentrate, except maybe for 

This question is addressed in an upcoming article [6]. 
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